Anisotropies of the lower and upper critical fields in MgB 2 single crystals 
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The temperature dependence of the London penetration depth (A) and coherence length (£) has 
been deduced from Hall probe magnetization measurements in high quality MgB2 single crystals in 
the two main crystallographic directions. We show that, in contrast to conventional superconductors, 
MgB2 is characterized by two different anisotropy parameters (T\ = X c /X a b and I\ = £ab/£c) which 
strongly differ at low temperature and merge at T c . These results are in very good agreement with 
recent calculations in weakly coupled two bands suprerconductors (Phys. Rev. B, 66, 020509(R) 
(2002). 
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The recent discovery of superconductivity in Mag- 
nesium Diboride [1] has been the starting point of a 
large number of theoretical and experimental investiga- 
tions. The coexistence of two superconducting energy 
gaps with different anisotropies [2] is expected to give rise 
to very peculiar physical properties in this system [3-5]. 
Among those, the strong temperature dependence of the 
anisotropy of the upper critical field r# c2 — H c i\\a,bl H-ci\\c 
is now well established [6,7] (H C 2\\ a b and H c2 \\ c being 
the upper critical fields for magnetic fields parallel to 
the two main crystallograpic directions, ab— planes and 
c— direction, respectively). 

"Classical" superconductors can be characterized by 
one single anisotropy parameter : F = £ a b/£c — A c /A a f, 
(where £ and A are the superconducting coherence and 
penetration depths in signed crystallographic orienta- 
tions) whatever the anisotropy of the superconducting 
gap. However, it has been suggested that £, a b/£,c could dif- 
fer considerably from A c / A a h at low temperature in MgB2 
due to the presence of two superconducting energy gaps 
with different anisotropies [4]. Even though recent neu- 
tron scattering data have confirmed that A c /A a h ~ 1 at 
low temperature is indeed very different from £ a b /£ c ~ 5 
[8], the temperature dependence of T\ = \ c /\ a b still 
had to be determined. In this letter, we confirm that 
P? = iab/£,c (= ^ Hc2 ) strongly differs from T\ at low 
temperature and show for the first time that and T\ 
present different temperature dependencies : whereas 
decreases from ~ 5 to ~ 2 (close to T c ), T\ increases with 
temperature tending from ~ 1.4 towards for T — > T c . 

Magnetization measurements have been performed on 
high quality single crystal [9] (T c w 36.5 K) showing 
flat surfaces of typical dimensions : 100 * 100 * 25/im 3 
using a Hall probe magnetometer. The main surface 
(i.e. the ab— planes) of the sample has been placed ei- 
ther parallel to or perpendicular to the surface of the 



Hall probes in order to measure the magnetization for 
H\\c and H\\ab respectively (see sketches in the inset of 
Fig. la). The alignement of the external field with the 
main crystallographic axis has been obtained by slightly 
rotating the ensemble in order to get the maximum and 
mimimum H C 2 values for H\\ab and H\\c respectively. A 
typical magnetization loop is presented in the inset of 
Figlb for H\\c and T — 20 K. As shown the irreversible 
part of the magnetization loops is rather small (corre- 
sponding to critical current values on the order of 10 3 
A/cm 2 at low T and low H) and the reversible mag- 
netization could thus easily be obtained assuming that 
M rev = (M up + M down )/2 where M up and M down are 
the magnetization for increasing and decreasing magnetic 
fields respectively (Bean critical state model). Typical 
curves are displayed in Fig.l for T > 26K. For lower 
temperatures, the upper critical field values for H\\ab be- 
came larger than our maximum field (3 T) but the lower 
critical field could still be deduced easily in both direc- 
tions down to the lowest temperatures (5 K). 

For H\\c i.e. perpendicular to the platelet, important 
demagnetizing effect come into play and the magnetic 
field has been rescaled to : H = Hq — N c M rev where 
Hq is the external field and N c the demagnetizing factor. 
iV c can be estimated assuming that the sample is an el- 
lipsoid of thickness t ~ 25/im and width w ~ lOO^tm to 
N c = 1 — irt/2w w 0.6. This value has been confirmed by 
the fact that, after correction for this demagnetization 
effect, the magnetization curves present a nearly vertical 
slope for H — H c \. As pointed out by Zeldov et al. [10] 
this value might however be overestimated by a factor 
on the order of y 'w/t due to the presence of geometrical 
barriers but it is important to note that, even though 
the absolute values of -ff c i|| c - and i? c i|| a {, - are directly 
related to the choice of the demagnetizing factor, their 
temperature dependence is absolutely not affected by this 
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choice. 
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FIG. 1. Magnetic field dependence of the reversible mag- 
netization for H\\c (a) and H\\ab (b) at various temperatures. 
In the inset : (a): schematical drawings of the geometry used 
for both directions, and (b): complete magnetization loop at 
T = 20 K for H\\c. 



The value of N ab has then been set to N ab « 1 - N c /2. 
In this direction, the effect of Bean Livingston barriers 
can also lead to a significant overestimation of H c i (up 
to a factor K/ln{n) where n is the Ginzburg -Landau 
parameter, the penetration field is then equal to the 
thermodynamic field instead of the lower critical field). 
Those effects play for instance a significant role in high T c 
cuprates due to their very high k values [11]. They result 
from the competiting effect between an attracting image 
force to the surface and a repulsive one arising from the 
interaction between the vortices and the shielding cur- 
rents. This latter force is proportional to the magnetiza- 
tion and thus vanishes asM^O for descending magnetic 
fields. When the field further decreases nothing prevents 
the vortices to leave the sample and M remains close to 
zero. This effect is highly sensitive to the quality of the 
surface and we did never observe this characteristic be- 
haviour in our magnetization loops. On the contary the 
irreversible part of the magnetization remained small in 
the entire temperature range (being for instance on the 
order of 15% at 30 K) and we thus assumed that M rev 
is equal to the average magnetization (i.e. that pinning 
mainly arises from bulk defects). 



The corresponding M rev versus H curves are displayed 
in Fig.l. Note that the y-axis has been slightly rescaled 
in order to get a —1 slope in the Meissner state due to 
the fact that the sample did not completely recover the 
surface of the Hall probe (especially for H\\ab). Assum- 
ing that the lower critical magnetic field is equal to the 
penetration field (i.e. neglecting durfacc pinning effects, 
see discussion above), H c i could be easily determined 
from the well defined minima in the M rev (H) curves 
(see Fig.l). The upper critical magnetic field H C 2 has 
been defined as the onset of the diamagnetic response at 
M rev (H c2 ) = 0. 
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FIG. 2. Temperature dependence of the lower critical field 
for H\\c and H\\ab. In the inset : Temperature dependence 
of the upper critical field for the same directions. The dotted 
line is an extrapolation below 26 K deduced from our previous 
magnetotransport measurements [5] 

The temperature dependencies of the corresponding 
critical fields are displayed in Fig. 2. The inset shows 
the temperature dependence of the upper critical mag- 
netic field for H \\ c (solid circles) and H\\ab (open circles) 
deduced from our magnetization measurements. Those 
values have been extrapolated down to low temperature 
for H \\ab using our previous high field magneto-transport 
data performed on a sample from the same batch [6] 
(with slightly higher H c2 values, dotted line in the inset 
of Fig. 2). As previously observed [6,7], the temperature 
dependence of H c i is almost linear for H\\c and presents 
a positive curvature for H\\ab. As shown, the situation 
is very different for the lower critical field which is al- 
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most linear for H\\ab and reveals a negative curvature for 
H\\c at high temperatures. This unusual negative cur- 
vature of H c i\\ c close to T c has been previously observed 
in thin films and polycrystalline samples [12] and can be 
explained by the two-band Ginzburg-Landau theory [13] 
(the isotropic Ginzburg-Landau theory only allows for a 
linear temperature dependence of H c i close to T c [14]). 
As shown below, this difference in curvatures will lead to 
an increase of T\ for increasing temperature. 



gle crystals by Caplin et al. [19] but lower than 
those obtained on thin films and polycrystalline sampcs 
[20,21,12]. However, A is expected to be very sensitive 
to interband scatering and it has been shown [22] that 

X ab tV (°) ~ 100 nm whereas A^ ea ™(0) « 40 nm (and 
~ 300 nm an d A^ ea ™(0) w 40 nm) suggest- 
ing that our single crystals are close to the clean limit as 
further confirmed by our T\ value « 1.4 (r^ eon (0) ~ 1 
and T^ rty (0) ~ 3). This ratio is also in good agreement 
with the one recently deduced from neutron scattering 
measurements [8]. 
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FIG. 3. Temperature dependence of the anisotropy param- 
eters Tj, T\ and r_sr el (see text for details). The dotted line is 
an extrapolation of T 4 at low temperature deduced from our 
previous magnetotransport measurements [5]. 

The penetration depth has been deduced from H c \ 
through : 

H cl = ^ /(27T\ 2 ab ).(T H /T x ).(ln(k) + 0.5) (1) 

w ith k = K b/jab-, T ff = T A for H\\c and k « 
\J Xab^c/£,ab£,c, = 1 for H\\ab. The values of the coher- 
ence lengths £ c (0) ~ 5 nm and £ a &(0) « 10 nm have been 
determined from H C 2{T) and we hence got A c (0) « 80 
nm and A a &(0) « 60 nm in good agreement with the val- 
ues deduced from the slope of the magnetization curve : 
1/A 2 = [8n/$o\d(lJ,oM rev )/dln(H). Note that this slope 
rapidly decreases with field up to some characteristic field 
H*(T) being on the order of 0.3 T at low temperature. 
This field can probably be associated with the closing of 
the small gap, H* ~ H* 2 in good agreement with our 
previous point-contact spectroscopy experiment [15] and 
with estimation of Bouquet et al. [16] from specific heat 
measurements and Eskildsen et al. [17] from tunneling 
data. Note that it has been suggested that for uniaxial 
superconductors k could be approximated to A a fc/A c for 
all field directions [18] which would lead to A c m X a b ~ 60 
nm. 

The \ a b(0) value determined in this experiment is 
much larger than the one previously obtained on sin- 
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FIG. 4. Temperature dependence of d(l/\ 2 )/dT show- 
ing a maximum (resp. minimum) in the ab— planes (resp. 
c— direction). In the inset : temperature dependence of 1/A 2 
for the two main crystallographic directions. 

The temperature dependence of the as-declucecl T\ val- 
ues is shown in Fig. 3 as solid circles. T^(T) (open circles) 
has been deduced from our H c2 values above 26 K and ex- 
trapolated towards lower temperature using the H C 2(T) 
extrapolation shown in the inset of Fig. 2. As shown, 
in striking contrast with conventional one band super- 
conductors the anisotropy parameters T\(T) and T^(T) 
are very different at low temperatures. Moreover, they 
reveal opposite temperature dependencies and Fx (T) in- 
creases tending towards T^(T) for T — > T c where both 
anisotropics merges at r^(T c ) = Y\{T C ) ~ 2. This un- 
usual behavior is also in striking contrast with the one 
previously obtained by Caplin et al. [19] who suggested 
that = Y\ ~ 2 in the entire temperature range but is 
in a very good agreement with Kogan's [4] calculations 
of the anisotropy parameters for a weakly-coupled two- 
bands superconductor. 

At low temperature the anisotropy of the upper crit- 



3 



ical field is mainly related to the anisotropy of the 
Fermi velocities over the quasi 2D a sheet : w 
V < v ab ><T / ~ v c >a ~ 6. then decreases with tem- 
perature as the influence of the small, nearly isotropic, 
gap increases due to thermal mixing of the two gaps [4,5] . 
On the other hand, the anisotropy of A is related to the 
anisotropy of the Fermi velocities over the whole Fermi 
surface which is expected to be on the order of 1.1 in 
reasonable agreement with our experimental data. As 
pointed out by Kogan et al. [4] the two anisotropics have 
to merge at T c as they are then determined by the same 
"mass tensor". We have also reported in Fig. 3 the tem- 
perature dependence of Th c1 = ffci||c/-^ci||o6- As shown 
Tn al and T\ present very similar temperature dependen- 
cies (and only differ by a numerical factor on the order 
of 1.4) showing this dependence does not depend on the 
choice of k in Eq.(l). 

Finally we discuss the temperature dependence of the 
superfluid density oc 1/A 2 . The inset of Fig. 4 dis- 
plays (A(0)/A(T)) 2 as a function of temperature for the 
two main crystallographic directions. The dotted lines 
schematically represents the temperature dependencies 
of the superfluid densities in the case of two indepen- 
dent BCS superconducting n (for which T c would be on 
the order of 20 K) and a bands. At low temperatures, 
the T— dependence of 1/A 2 . and 1/A 2 fc are both deter- 
mined by the the band with the smallest gap (i.e. the 
7T— band) and T\ is only very weakly temperature de- 
pendent. As the temperature increases, the influence 
of the a— band sets in and an inflexion point is then 
clearly visible around 20K for l/A 2 fc (T) (see maximum 
for d{l/X 2 ab )/dT in Fig. 4) followed by a downward cur- 
vature up to T c . As pointed out in [22], this inflexion 
point is a direct consequence of the "superposition" of 
two BCS bands. The sharp kink that would appear for 
two independent bands is here smoothed out by inter- 
band coupling but still shows up as a clear maximum in 
the first derivative close to the "critical temperature" of 
the 7T— band. A very similar behaviour has been previ- 
ously reported by Carrington et al. [21] for the superfluid 
density deduced from microwave measurements. 

It is important to note that this change in the curva- 
ture at T = 20 K is also visible for 1/A 2 (see mimimum 
for d{\/\ 2 c )/dT in Fig. 4). In this direction, the pene- 
tration depth is mainly determined by the 7r— band up 
to T c (in the clean limit) due to the very small value of 
the superfluid density of the a— band (i.e. strongly 2D 
character of the er— band). Despite the small value of the 
corresponding gap, superconductivity remains induced in 
the 7r band by intcrband coupling giving rise to the " tail" 
of 1/A 2 at high temperature which shows up as a lower 
<i(l/A 2 )/<iT value. Those temperature dependencies are 
qualitatively in good agreement with the recent calcula- 
tions by Golubov et al. [22]. However, the details of those 
dependencies strongly depend on intraband scattering as 



well as on the relative contributions of the two bands. 

In conclusion, we have given direct experimental ev- 
idence that the anisotropy of MgB 2 is characterized by 
two different parameters (T\ and T^) which differs not 
only in absolute values but also in their temperature de- 
pendencies in striking contrast to conventional type II 
superconductors. The as-deduced temperature depen- 
dencies of those two anisotropy parameters are in perfect 
agreement with theoretical predictions for weakly cou- 
pled two bands superconductor [4] . The presence of two 
different bands is enforced by the inflexion point observed 
around 20 K in the temperature dependence of 1/A 2 in 
both crystallographic directions. 

The work of P.Sz. was partially supported by the Slo- 
vak Science and Technology Assistance Agency under 
contract No.APVT-51-020102. This work in Korea was 
supported by the Ministry of Science and Technology of 
Korea through the Creative Research Initiative Program. 



[1] J. Nagamatsu et al., Nature 410,63 (2001). 

[2] P. Szabo et al., Phys. Rev. Lett. 87, 137005 (2001), 

F.Giubileo et al. Phys. Rev. Lett. 87, 177008 (2001), 

F. Bouquet et al. Europhys. Lett. 56, 856 (2001). 
[3] A. Y. Liu et al., Phys. Rev. Lett. 87, 087005 (2001). 
[4] V.G.Kogan, Phys. Rev. B 66 020509(R) (2002); P. Mi- 

ranovic et al. J. Phys. Soc. Jpn. 72 221 (2003). 
[5] T.Dahm and N.Schopol, Phys. Rev. Lett. 91 017001 

(2003). 

[6] L. Lyard et al. Phys. Rev. B. 66, 180502(R) (2002). 
[7] S.L. Bud'ko et al. Phys. Rev. B 64, 180506 (2001), 
M.Angst et al. Phys. Rev. Lett. 88, 167004 (2002), 
U.Welp et al. Phys. Rev. B 67, 012505 (2003). 
[8] R.Cubitt et al. Phys. Rev. Lett. 90, 157002 (2003). 
[9] K. H.P.Kim et al. Phys. Rev. B 65, 100510(R) (2002). 
[10] E.Zeldov et al. Phys. Rev. Lett. 73, 1428 (1994). 
[11] M.konczykowski et al. Phys. Rev. B 43 , 13707 (1991). 
[12] For a review see : C. Buzea et al. Supercond. Sci. Tech- 

nol. 14 R115 (2001) 
[13] I.N.Askerzade et al. Supercond. Sci. Techno. 15 L17 
(2002) 

[14] A.A.Abrikosov, in book Fundamentals of the theory of 

metals, North-Holland, Amsterdam (1988). 
[15] PSamuely et al. Physica C 385, 154 (2003). 
[16] F. Bouquet et al. Phys. Rev. Lett. 89, 257001 (2002). 
[17] M.R.Eskildsen et al. Phys. Rev. Lett. 89 187003 (2002). 
[18] A.B. Balatskii et al. Sov. Phys. JETP 63 866 (1986). 
[19] A.D.Caplin et al. Supercond. Sci. Technol. 16 173 (2003). 
[20] X.H.Chen et al. Phys. Rev. B 64, 172501 (2001), 

C.Niedermayer et al. Phys. Rev. B 65 094512 (2002). 
[21] F.Manzano et al. Phys. Rev. Lett. 88 047002 (2002), 

A. Carrington and F.Manzano, Physica C 385, 205 

(2003). 

[22] A. A. Golubov at al. Phys. Rev. B. 66, 054524 (2002). 



4 



